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Abstract. We propose a simple picture for the occurrence of superconductivity 
and the pressure dependence of the superconducting critical temperature, TgQ, 
in ZrZn2. According to our hypothesis the pairing potential is independent of 
pressure, but the exchange splitting, Exc, leads to a pressure dependence in the 
(spin dependent) density of states (DOS) at the Fermi level, D^iep). Assuming 
p-wave pairing Tgc is dependent on Daiep) which ensures that, in the absence 
of non-magnetic impurities, Tgc decreases as pressure is applied until it reaches 
a minimum in the paramagnetic state. Disorder reduces this minimum to zero, 
this gives the illusion that the superconductivity disappears at the same pressure 
as ferromagnetism does. 
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The coexistence of ferromagnetism and superconductivity is a problem of long 
standing and general interest. Thus its recent discovery in UGe2 1 , ZrZn2 [2] 
and URhGe j3j is attracting considerable attention. In particular its occurrence 
in ZrZn2 is intriguing because, at ambient pressure, ZrZn2 is a weak ferromagnet 
{TpM ~ 28.5 K) and by the application of pressure it can be tuned through a quantum 
critical point (QCP) [Pc ~ 21 kbar) to become a paramagnetic metal 0. This revives 
an old suggestion of Fay and Appel 0]. These authors calculated Tsc mediated by 
paramagnons in a McMillan like formalism and found that there is superconductivity 
in the (triplet) Ai channel on both sides of the QCP. However, while the broad 
description of Fay and Appel agrees with the observations the fine details do not. In 
Fay and Appel's theory the transition temperature goes to zero at the QCP, Tsc then 
rises to a (local) maxima as the model is tuned away from criticality (experimentally 
this corresponds to pressure being varied away from Pc)- Tsc then falls again away 
from the QCP. They also predicted that Tsc would be approximately the same 
magnitude in both the ferromagnetic and paramagnetic sides of the phase diagram 
(although slightly higher on the paramagnetic side). 

When superconductivity was observed in ZrZn2 it was only seen on the 
ferromagnetic side of the phase diagram and further the maximum in Tsc was observed 
at ambient pressure The experiments show a monotonic decrease in Tsc with 
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pressure until about 18 kbar. No data has been published in the pressure range 18 
kbar < P < 22 kbar, therefore one cannot ascertain where in this range Tsc falls to 
zero. Several groups have attempted to explain this either by revisiting the theory of 
Fay and Appel and examining specific coupling mechanisms |S] or by considering the 
problem within a Ginzburg-Landau formalism [7] . Both of these groups predicted that 
Tsc goes to zero at Pc- Here we will present a simple alternative to these scenarios 
for the coexistence of superconductivity and ferromagnetism in ZrZn2. We will show 
that no variation in the coupling constant with pressure (i.e. proximity to the QCP) 
is required to explain the experiments due to the natural enhancement of the Ai 
phase transition temperature in the ferromagnetic phase and the extreme sensitivity 
of the triplet pairing states to scattering from non- magnetic impurities. On the other 
hand our arguments rely on ZrZn2 being a rare example of a Stoner ferromagnet. 
Interestingly in our consideration the QCP plays no special role and the pressure at 
which superconductivity disappears is predicted to be strongly sample dependent. 

We wish to consider the problem via the simplest model which has the possibility 
of illustrating the relevant physical phenomena: triplet superconductivity and 
ferromagnetism. To this end we study a one band Hubbard model with an effective, 
attractive, pairwise, nearest neighbour interaction, Uija-a-' ■ We allow ferromagnetism 
to enter via the Stoner model which appears to be in good agreement with the observed 
behaviour of the ferromagnetic phase of ZrZn2 |2 |H1 IHl Thus we study the 

consequences of the following Hamiltonian: 

Ti. — — ^ij^\(y^j(y + 2 5Z ^yo'0''^lcr^*'T^]o-'^Jo-' ~ CrExccl^Cja- (1) 
ija ijaa' icr 

where c'^J are the usual annihilation (creation) operators for electrons with spin 
(7 = ±1 occupying a tight binding orbital centred on the lattice site labelled by i. 
To render the model tractable we assume that the sites i form a simple cubic lattice. 

For a general triplet pairing state in a field gap equations cannot be derived in 
the same way as they can for the zero field case ^2 However, if we specialise 
to the case of equal spin pairing (ESP) and neglect the action of the dipolar field 
on the orbital motion of the electrons a remarkable simplification occurs as shown in 
reference [T^I the gap equations are 

^-(^) = - E - 2/-.v) (2) 

k' 

where the quasiparticle spectrum is given by 

Eka = v'(efc-A*-fTi^.c)' + |A,^(fc)|2. (3) 

Equations and © have several surprising features. Firstly, there is a complete 
decoupling of the two spin states; even in the presence of Cooper pairing. Secondly, 
the exchange splitting enters only in the role of a chemical potential, but with opposite 
signs for the two spin states. It must be stressed that these results are only valid for 
ESP states. However, the large exchange splitting in a ferromagnet probably precludes 
opposite spin pairing (OSP) states (that is singlet states, = Q triplet states or 
even the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, which is only stable for 
moderate exchange splittings). In short (0) and lO are just the Hartree-Fock-Gorkov 
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Figure 1. The DOS from the LDA calculations for ZrZn2 ( ) by Santi and 

coworkers 021 EH and our tight binding fit to the LDA DOS ( ). Inset: The 

experimental data llb l is consistent with a fit to TgQ{P = 0) ~ 1.2 K from the 
Abrikosov-Gorkov formula l5t . 



approximations restricted to the ESP sector of the theory. For Uija-a' = UijS^^i no 
such restriction is required as OSP states are not a possibihty. 

As T ^ Tsci |Acrcr'(fc)| and we find the hnearised gap equations: 

A„„{k) = > — ^ tanh — — A^^ A; . 4 

2[e{k') - fj, - (jExc) \ ^ksT J 

Before solving these equations numerically we must first choose the hoping integrals, 
tij, and the coupling constants Uijaa'- We fit the hoping integrals with an on site 
{tii = fj.) and nearest neighbour terms only so as to give the same relative density of 
states in the region of the Fermi level as is found in ab initio band structure calculations 
The DOS from our fit is compared with that found in the ab initio calculation 
in figure n Evidently, our one band model cannot reproduce the complex behaviour 
of the DOS over the full energy range of the Zr related d-band, nevertheless, within an 
energy range of 20 meV about the Fermi energy it does do so. Since Uijaa' depends 
on only one parameter: Uij = U for sites i and j being nearest neighbours (Uij = 
otherwise) it can be determined by reference to the measured Tgc for clean ZrZn2. 
This is hampered by the extreme sensitivity of triplet pairing to scattering from non- 
magnetic impurities |14[ [T5] and by the lack of data. Using experimental estimates of 
the residual resistivity, p, we find that what little data there is [IB] is consistent with a 
clean superconducting critical temperature, at ambient pressure, Tg(j{P = 0) ^ 1.2 K 
as shown in the inset to figure ^ 

We solved numerically with U = 0.88i on a k-space integration mesh of 10^ 
points. Such a fine integration mesh is required to accurately reproduce the DOS. Our 
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method (implicitly) requires an accurate calculation of the D„{ep) as we are varying 
the exchange splitting and thus we are changing D^lep), so any errors in evaluating 
Daisp) will lead to significant errors in our calculation of the variation of Tsc with 

The results of our numerical calculations are shown in figure El The Ai 
phase displays superconductivity in only the tt channel. In the Ai phase d{k) ~ 
(fc^ + ifcj,, i(fc^ + ifcj^), 0), where d{k) is the usual BW vector order parameter for 
triplet superconductivity 0. The A2 phase corresponds to superconductivity in both 
ESP channels but with different amplitudes in the two channels so that d{k) ~ 
(kx + iky , iK{k,j; +ifcj^),0) where < k < 1. The A phase, which is only stable in 
zero exchange spitting, has the same pairing amplitude for both of the ESP states, 
corresponding to the order parameter d{k) ~ {k^ + iky, 0, 0). Because of the complete 
separation of the up and down sheets the linearised gap equations Q can be used 
to calculate the lower transition temperature, T^q. This transition represents the 
formation of a superconducting state in the minority spin state. Our phase diagram 
of course assumes that no other phase transitions occur. The existence of the Ai and 
A2 phases only is consistent with a general symmetry analysis \IJ\ and the Ginzburg- 
Landau expansion of the free energy |12[ [TH| . 

To make contact with experiment we must allow for the strong dependence of the 
superconducting transition temperature of a p-wave superconductor on non-magnetic 
impurity scattering |14[ IT^ . This is done via the Abrikosov-Gorkov formula: 

ln(^)^J'- + ^ )-Jl) (5) 

\Tsc J ^ ^^TtrkBTsc ) ^V2y ^' 

where Ttr is the quasiparticle lifetime as measured in transport experiments and '4){x) 
is the digamma function. Note that we do not need to worry about the Baltensperger- 
Sarma equation |ITj)ir2nj . which accounts for the reduction in the critical temperature of 
a superconductor due to exchange splitting, as this is only valid for OSP states. Thus, 
the Abrikosov-Gorkov formula can be used to calculate Tsc as a function of rt^, or 
equivalently p via the Drude formula. To make the most of the available experimental 
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Figure 3. The critical temperature of a ESP p-wave superconductor and tlie 
temperature of the A1-A2 transition as a function of exchange splitting in the 
presence of disorder. The curves correspond (from the top down) to ptr = 0, 
0.1 fiQcm, 0.2 /iQcm, 0.3 /iQcm, 0.4 fiQcm, 0.5 /jficm, 0.6 fiQcm, 0.7 fiQcm, 
0.8 jiQcm and 0.9 jiflcm. 



data we used the Abrikosov-Gorkov formula ^ to investigate the effect of disorder 
on the above phase diagram (see figure OJ. 

The final step needed to make a direct comparison with measurements of Tsc 
as a function of pressure is to note that, experimentally, the Curie temperature is a 
Hnear function of pressure Namely, Tfm{P) = Tfm{0) (1 - P/Pc)- The zero 

temperature magnetisation is also linear in pressure |S] and thus proportional to TpM 
(as is predicted by the Stoner model), giving M(P,T = 0) = M(0,0) (l-P/Pc)- 
For a Stoner ferromagnet, the magnetisation is linearly dependent on the exchange 
splitting and hence 



We now invoke the fact that TpM S> Tsc which implies that E^dP, T — Tsc) ~ 
Exc{P,T = 0). Thus we can map the results of Tsc{Exc) (shown in figure |3Jl onto 
Tsc{P) which we show in figure^l It can be seen that although quantitative agreement 
with experiment is not achieved, the general features of experiment are reproduced by 
several of the curves. 

Thus we conclude that we have demonstrated the viability of the following simple 
picture. Irrespective of the mechanism of pairing and exchange splitting ZrZn2 is 
a p-wave superconductor with a low Tsc ~ 1.2 K (at ambient pressure). This 
superconductivity is not observed in the paramagnetic phase of currently available 
samples due to disorder. However, the exchange field enhances Tsc of an Ai p-wave 
state and this is the cause of the observed superconductivity in the ferromagnetic 
phase. Experiment suggests PIIHIEIEI ZrZn2 is a rare Stoner ferromagnet for which 
the exchange splitting is proportional to the magnetic order parameter. Thus when 
P > Pc and therefore TpM = M = Exc = there is no measurable superconductivity. 
However, improvement in sample quality will lead to a lowering of the residual 




P<Pc 
P>Pc- 



(6) 



Letter to the Editor 



6 




Figure 4. The critical temperature of a ESP p-wave superconductor and the 
temperature of the A1-A2 transition as a function of pressure in the presence of 
disorder. The theoretical curves are scaled so that Tsc = 1 at ambient pressure. 
The experimental data, taken from was scaled in the same way after straight 
line had been fitted to the data. The curves correspond (moving from top right 
to bottom left) to ptr = 0, 0.1 fiQcm, 0.2 fiQcm, 0.3 fiQcm, 0.4 fiQcm, 0.5 fiQcm, 
0.6 ixQcm, 0.7 /if2cm, 0.8 fiQcm and 0.9 fiQcm. 



resistivity and thus presents the possibility of the observation of superconductivity 
in the paramagnetic state (as is demonstrated by the curves with ptr < 0.3 ^Ocm in 
figure 2J). This explanation is consistent with and lends microscopic support to the 
more phenomenological arguments of Walker and Samokhin [7] and Mineev J7] . 
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